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EFFECTS OF ASYMMETRIC FLOWS IN SOLAR CONVECTION ON OSCILLATION MODES 
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ABSTRACT 

Many helioseismic measurements suffer from substantial systematic errors. A particularly frustrat- 
ing one is that time-distance measurements suffer from a large center to limb effect which looks very 
similar to the finite light tra vel time, except t hat the magnitude depends on the observable used and 
can have the opposite sign (jZhao et al.ll2012[ ). Th is has frust r ated a ttempts to determine the deep 
meridional flow in the solar convection zone, with I Zhao et "all ()2012D applying an ad hoc correction 
with little physical basis to correct the data. In this letter we propose that part of this effect can be 
explained by the highly asymmetrical nature of the solar granulation which results in what appears to 
the oscillation modes as a net radial flow, thereby imparting a phase shift on the modes as a function 
of observing height and thus heliocentric angle. 

Subject headings: Sun: helioseismology — Sun: oscillations — Sun: granulation 
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1. INTRODUCTION 

The meridional flow in the solar convection zone plays 
a key role in many solar dynamo models and an accu- 
rate measurement of the flow with depth and latitude 
would thus be invaluable for constraining solar dynamo 
models. There have been numerous attempts to obtain 
estimates using a variety of techniques , including time- 
distance (jGilesI 120001: iZhao et al.! 120121. in th e following 
Z12) , ring diagrams (ISchou & Bogart 1998; Haber et al.l 
I2002f ). normal modes (|Woodard et al.ii2012i . in the fol- 
lowin g W12), and supergranulation studies (jHathawavl 
1201 If ). While these studies have given reasonable num- 
bers near the surface, they have suffered from large and 
unexplained systematic errors, preventing us from ob- 
taining reliable numbers throughout the convection zone. 

A meridiona l counter-ce ll at high latitude was 
first found _by Haber et al.l (120021 ) : subsequen t work 



(| Gonzalez Hern andez et al.l 120061 iZaatri et all I2006D 

found this to be a periodic phenomena t ightly corre- 
lated with the solar inclination angle (i?o)- IZaatri et al.l 
|2P06), concluding that the counter-cells were likely 
spurious, app li ed a correction to remove them. 
IBraun fc BirchI ()2008n found that North-South travel 
times differed depending on the heliocentric longitude 
at which the measurement was performed. Most re- 
cently, Z12 used data from the Helioseismic and Mag- 
netic Imager (HMI) to measure East- West travel time 
shifts along the equator and North-South travel time 
shifts along the central meridian in four different observ- 
ables: continuum intensity, line core intensity, line depth 
(continuum minus line core) , and Doppler velocity. They 
found large E-W travel time shifts along the equator, and 
found that they were quite different in different observ- 
ables. 

Although Z12 did not provide an explanation for the 
source of the error, they treated it as a heliocentric angle 
dependent phase or time shift. Using this assumption, 
they used the E-W travel time anomalies to correct the 
N-S travel time shifts — this brought the four different 
observables into good agreement, which was encouraging. 
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Similarly, it was noted by W12 that a radially varying 
phase of the eigenfunctions coupled with the variation in 
the height of formation with changing heliocentric angle 
might be an explanation, but again no source of such a 
phase variation was identified. 

It is evident that helioseismic observations should suf- 
fer fr om a phase error as a function of heliocentric angle 

— as iDuvall fc Hanasoei ([1009) pointed out, the light 
travel time for an observation at disk center is differ- 
ent than for an observation near the limb by roughly 
2s. When analyzing travel time residuals, however, they 
found the correct magnitude but opposite sign. They 
conlcuded that the measured travel times suffer from a 
systematic effect with twice the magnitude and opposite 
sign as the expected light travel t i me eff ect. Similar num- 
bers were found by ISchou et al.l ()2012f ). who found that 
a travel time error of 2-3s at a heliocentric angle of 60° 
could explain their results. 

In general, standing acoustic waves should have a con- 
stant phase with height. Many effects that are known to 
be poorly modeled or neglected entirely do not change 
this. What is required to add a phase shift with height is 
an asymmetric effect, e.g. an effect which knows whether 
a wave is traveling upwards or downwards. In the Sun, 
of course, the modes we observe are not purely standing, 
and mode damping due to, for example, non-adiabatic ef- 
fects will have an effect. In this work we are considering 
low frequencies, however, so we neglect this. 

Here we suggest that a phase variation arises from the 
large asymmetry between the upflows and downflows in 
the convection near the solar surface. In particular, the 
broad upflows and narrow downflows give rise to net ver- 
tical flows when horizontally averaged over length scales 
much smaller than the aco ustic modes. As was shown by 
I Cough fc HindmanI ()2010[ ) in the context of meridional 
flows, a flow introduces a phase shift in acoustic modes 

— one would expect the vertical convective flows in the 
Sun to introduce a phase shift with height in the solar 
atmosphere. 

In the following discussion, we extract complex eigen- 
functions from a detailed numerical simulation of convec- 
tion to compute the phase delay as a function of height. 
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Figure 1. Mach number as a function of depth from a model by 
I Stein et ai.l 120091 ). The area average of the flow divided by the 
area average of the sound speed is used. To exclude the buffer 
zone used in the model the outer 5 grid points were excluded. The 
Mach number below IMm is negligible. 

We attempt to explain this phase shift as being due to 
the spatially averaged vertical flows by computing phases 
shifts to theoretical eigenfunctions with an imposed ver- 
tical flow and comparing these to the numerical eigen- 
functions from the simulation. We use the phase shifts 
from the simulation data to predict systematic effects in 
travel-time measurements as a function of distance from 
disk center, and compare these effects to the observed 
discrepancies in solar data. Finally, we discuss some of 
the shortcomings of our models, how they might be im- 
proved, and discuss other observable consequences. 

2. MODELS 

We use a non-rn agnetic simulation of solar convection 
(jStein et al.ll2Q09[ )PT Convection in the outer layers of 
the solar envelope and atmosphere are simulated in a 
small Cartesian box. We show the horizontally averaged 
vertical Mach number in Figure [TJ In this work we as- 
sume that the modes 'see' only this horizontal average. 
For low degree modes, this is likely to be an adequate 
approximation for our purposes. 

2.1. Simulation Results 

The phase shifts of the eigenfunctions can be estimated 
by extracting them from a convection simulation, as pre- 
viously done (jStein fc Nordlundll2001[) . but this time in- 
cluding the imaginary component. This is easily done 
by horizontally averaging the vertical velocity (to isolate 
radial modes) from such a simulation and taking a tem- 
poral Fourier transform at each depth. The phase shifts 
can be converted to time delays by dividing by the an- 
gular frequency uj. Results of such an analysis are shown 
in Figure [2] for the two lowest frequency modes (those at 
higher n are less well defined). 

2.2. Explaining the Simulation Results: Theoretical 
Eigenfunctions 

What is the cause of the phase shift we see in the stand- 
ing waves in the simulation box? To get a crude estimate 
of the effect of the vertical fl ows, we start by con s iderin g 
the wind-in-a-pipe model of lGough fc Hindmanl (|2010( ). 

^ The model can be found at |http: //sha. Stanford. edu/stein_sim/| 



200 




-50 I \ \ \ I \ \ 

-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
Depth (Mm) 

Figure 2. Time delay as a function of depth from the simulation. 
Solid line is for lowest n at 1142/xHz, dashed for the second lowest 
n at 1761^Hz. 

Section 3. Manipulating their equations and extending 
their Equation (9) to allow all variables to depend on 
position shows that the time shift introduced by a flow 
is given by 

Ar(a;) = / ^^dx = / ^^dx, 1 

Jxl C{XY C(x) 

where x is the distance from the end of the pipe. The re- 
sults of this calculation are shown in Figure [H It will be 
noted that, while some features of the numerical eigen- 
functions are qualitatively reproduced (i.e. a negative 
slope in the atmosphere and a steeper positive slope just 
below the photosphere) , there is a quantitative disagree- 
ment of more than an order of magnitude. This is not 
surprising, however, as this toy model assumes that the 
background model varies slowly compared to the wave- 
length, and that is manifestly not the case here. 

A somewhat more sophisticated approach is to com- 
pute the oscillation eigenfunctions for a stellar envelope 
with a specified vertical flow U{r). To do this, we per- 
turb the equations of continuity, motion, an d energy in 
the usual way (following, e.g., lUnno et all [T989il . but 
with the added velocity term in the base equations. For 
simplicity, we consider only radial modes in this work, 
but the generalization to non-radial modes is straight- 
forward. When a vertical flow is present, the solutions 
become complex, and we can compute a phase delay as 
a function of height. 

By including only a horizontally averaged vertical fiow, 
we are of course involving a certain physical inconsistency 
— that is, a horizontally invariant vertical flow is not con- 
sistent with a one dimensional background model, which 
requires a zero net mass fiux. Because we are computing 
eigenfunctions in the linear perturbative regime, this is 
not necessarily an unreasonable inconsistency to accept, 
but it does require that we employ certain assumptions. 
As noted above, we assume here that the modes 'see' the 
horizontal spatial average of the convective fiows. One 
consequence of this assumption is that we assume that 
variations on convective length scales (horizontally) do 
not affect the acoustic modes. For radial modes this is 
reasonable. It also also assumes that the correlations be- 
tween flows and thermodynamic quantities — say, den- 
sity — do not have an effect on the phases of these modes. 
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Figure 3. Predicted effect of vertical velocity on mode eigenfunc- 
tions. The dotted line shows the result of integrating Equation [l] 
multiplied by a factor of 10 to that it can be seen on this scale. The 
solid and dashed lines show the predictions from a more detailed 
calculation, at 1142 fiHz and 1761 /xHz, respectively. 

Using the horizontally averaged thermodynamic quan- 
tities and vertical flow from the simulation box, we in- 
tegrate the oscillation equations for radial modes with 
frequencies of 1142/iHz, and 1761/iHz, which may be 
directly compared with the eigenfunction phases shifts 
shown in Figure [51 The calculated eigenfunction phase 
shifts are shown in Figure [3l As can be seen, there is 
a general qualitative agreement between the two sets of 
phase delays. Most prominently, we see a strong positive 
phase shift just below the photosphere, though we calcu- 
late a somewhat smaller shift, and also find the peak to 
be much broader below the surface. In the atmosphere 
itself we find a fairly constant phase, or a slight phase 
lag with height. This is not entirely in agreement with 
the numerical eigenfunctions, but there is some uncer- 
tainty in the numerical eigenfunctions due to noise. We 
consider the agreement to be sufficiently good to con- 
clude that the phase shifts we observe in the convection 
simulations are due to the vertical convective flows. 

2.3. Effects on observations 

The HMI instrument measures, among other things, 
continuum intensity, line core intensity, and line-of-sight 
velocity. These observables are produced at a ran ge of 
heights in the solar atmosphere (jFleck et al.ll20ljT ). and 
can be represented as the convolution of some contribu- 
tion function of height and the actual variation of the 
quantity in the solar atmosphere. A contribution func- 
tion for t he HMI line core in tensity measurements can be 
found in lFleck et al.l (|2011l Figure 2); the other observ- 
ables have contribution functions that peak at different 
heights. The contribution function can be approximated 
by a Gaussian with a width of 250km. We will use this 
in the work that follows. A detailed study of the con- 
tribution function of the various HMI observables as a 
function of viewing angle is far beyond the scope of this 
letter. 

To relate the change with height to the change with 
viewing angle we need to determine the relationship be- 
tween these two quantities. Assuming that the atmo- 
sphere is isothermal and that the height of formation 
corresponds to the place where a certain column density 
of matter has been traversed from infinity, it can easily 
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Figure 4. Phase shifts as a function of viewing angle for different 
disk center formation heights, integrated over a Gaussian contri- 
bution function. The width of the Gaussi an is 250km, chosen to 
match the contribution function from Fleck et al.l I I2011I) . The top 
panel shows the phases at different viewing angles; the bottom 
panel shows the derivative of the phases with viewing angle. 

be shown that the change in height with angle is given by 
Ah = — H \og{cos{9)) , where H «120 km is the density 
scale height and 9 is the angle between vertical and the 
line of sight. At 6* = 60° this corresponds to about 80km. 

In Figure [4] we show the time delays integrated over 
a Gaussian contribution function as a function of view- 
ing angle for various disk-center formation heights. We 
also show the derivatives with respect to angular distance 
from disk center of the time delays. These derivatives 
give the travel time differences one would measure in the 
limit of infinitely small apperture size, and can provide a 
good approximation for the apperture sizes used in Z12. 
For direct comparison with Z12, the values of our deriva- 
tives must be multiplied by the apperture size used in 
the measurement. 

3. DISCUSSION 

We have explored the effects of vertical convective flows 
on solar acoustic modes. Can they explain the system- 
atics that have been observed in solar data, for exam- 
ple the effect found in Z12? In that work, the authors 
showed East- West travel time differences across the the 
solar equator for different HMI observables (line-of-sight 
velocity, continuum intensity, line core intensity, and line 
depth). They found that all observables exhibited sys- 
tematic effects along the equator that they believed to 
be spurious, but that each different observables had dif- 
ferent signatures. Comparing their Figure 2 with the 
bottom panel of our Figure |4l we note some striking sim- 
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ilarities. First, the continuum shows the largest anomaly 
— as does the 0km height measurement we predict, and 
with the same sign (negative). The Doppler signal, with 
a peak in the contribution function at approximately 
150km, shows a much smaller anomaly but still the same 
sign. We predict the same thing (note that our 150km 
signal changes sign, but only above 60° from disk cen- 
ter, which is further than the Z12 results go. Finally 
the line core intensity, with the highest formation height 
(250km), shows an anomaly with the opposite sign, and 
this is matched by our results. Furthermore, the Z12 
anomalies show turnovers as the distance from disk cen- 
ter gets large, and we show the same. 

The quantitative agreement is not as promising, how- 
ever. In particular, the magnitude of the effect we pre- 
dict, by multiplying the derivative of the phase shifts 
by the aperture size used in the measurements, is too 
small — at most a few seconds, as opposed to the more 
than 10 seconds Z12 find in continuum intensity. Fiiv- 
thermore, while we do predict the turnover observed, we 
do not accurately predict where this should happen. 

Several approximations may be identified, some of 
which are easily amenable to improvement and others 
not. First of all the use of the area weighted vertical ve- 
locity is unlikely to be correct. In reality the propagation 
of waves in a medium with small scale (^ wavelength) 
variations is a complicated issue and the velocities would 
likely have to be weighted in some other way. In particu- 
lar, variations in different quantities (significantly density 
and velocity) are correlated, which will need to be taken 
into account. 

Another problem is that the granulation is not static on 
the timescale of the oscillations and that there are likely 
to be interactions between the two. This should be well 
captured by the simulations, but is probably difficult to 
model accurately — at least for the present authors. 

Finally we have ignored complex issues of radiative 
transfer. In approximating the contribution function as 
a Gaussian, we have simplified the problem but done vi- 
olence to the actual physics. For the present purposes 
we consider it sufficient, but the exact shape of the ac- 
tual contribution functions (which are different both for 
different observables and for the same observable at dif- 
ferent viewing angles) will have significant quantitative 
effects. It is not likely, however, that the qualitative re- 
sults would be affected. At least this problem is well 
understood and has been addressed in detail for other 
purposes. 

In addition to the phase shift with observing angle, 
there are other possible observational consequences. Per- 
haps most obviously there will be a phase shift between 
different observables, such as continuum intensity and 
Doppler shift, which might be misinterpreted as a prop- 
agation or non-adiabaticity effect. 

Another effect is that even for the same observable 
there should be a phase difference between observation 
heights. This includes, e.g., Doppler shifts derived from 
different parts of the spectral line and even observations 
in the middle of granulation versus intergranular lanes. 

In addition to the observational effects, the fact that 
eigenfunctions can be determined with significant accu- 
racy from numerical simulations presents many opportu- 
nities, in particular ones involving the ability to deter- 
mine if the depth variation of both the real and imagi- 



nary parts match models. Potentially this could be used 
to test models of the interaction of waves with granu- 
lation and non-adiabatic effects, hopefully leading to a 
better understanding of the physics behind such things 
as the surface terms currently being applied in an ad hoc 
way in structure inversions. 

4. CONCLUSION 

We have shown that the effect of the vertical flows from 
convection in the outer solar convection zone and atmo- 
sphere do affect the quantities we observe in helioseismol- 
ogy. We have further shown that the systematic errors 
that have been observed can be qualitatively explained 
by this effect. We conclude that the ad hoc correction 
applied by Z12 is likely justified. 

A full quantitative prediction of the effects of the ver- 
tical flows requires a more sophisticated effort than that 
employed in this work. In particular, the modeling of the 
structure of the atmosphere must be very accurate and 
a proper treatment of the actual measurements we take 
must be done. This work, while non-trivial, is certainly 
feasible, and would be useful in addressing a number of 
different outstanding problems in helioseismology. 
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